UNUSUAL SYMMETRIES IN THE KUGEL-KHOMSKII HAMILTONIAN 



A. B. Harris 

Department of Physics and Astronomy, University of Pennsylvania, Philadelphia, PA 19104 



o 
o 



(N 



h 

CO 

c3 



I 

c 

o 
o 



> 

(N 
m 
o 
m 
O 

03 



o 
o 



X 



Taner Yildirim 

NIST Center for Neutron Research, National Institute of Standards and Technology, Gaithersburg, MD 20899 

Amnon Aharony, Ora Entin-Wohlman, and I. Ya. Korenblit 

School of Physics and Astronomy, Raymond and Beverly Sackler 
Faculty of Exact Sciences, Tel Aviv University, Tel Aviv 69978, Israel 
(Dated: February 2, 2008) 

The Kugel-Khomskii Hamiltonian for cubic titanates describes spin and orbital superexchange 
interactions between d 1 ions having three-fold degenerate t^ g orbitals. Since orbitals do not couple 
along "inactive" axes, perpendicular to the orbital planes, the total number of electrons in \a) 
orbitals in any such plane and the corresponding total spin are both conserved. A Mermin- Wagner 
construction shows that there is no long-range spin ordering at nonzero temperatures. Inclusion of 
spin-orbit coupling allows such ordering, but even then the excitation spectrum is gapless due to a 
continuous symmetry. Thus, the observed order and gap require more symmetry breaking terms. 

PACS numbers: 75.10.-b,71.27.+a,75.30.Et,75.30.Gw 



High temperature superconductivity^] and colossal 
magnetoresistanceQ sparked much recent interest in the 
magnetic properties of transition metal oxides, particu- 
larly those with orbital degeneracy Q. In many tran- 
sition metal oxides, the d electrons are localized due to 
the large on-site Coulomb interaction, U. Assuming a 
simple Hubbard model, with a typical nearest neighbor 
(nn) hopping energy t, the low energy behavior can be 
described by an effective superexchange model, which in- 
volves only nn spin and orbital coupling, with energies of 
order e = t 2 /U. In cubic oxide perovskites, the crys- 
tal field of the surrounding oxygen octahedra splits the 
d-orbitals into a two-fold degenerate e g and a three-fold 
degenerate t^g manifold. In most cases, these degenera- 
cies arc further lifted by a cooperative Jahn- Teller (JT) 
distortion^, and the low energy physics is well described 
by an effective superexchange spin-only model [E IE • 
However, some cubic perovskites, such as the titanates 
(RTi03, where R= La, Y, etc), have only a small JT 
distortion^, in spite of the orbital degeneracy 9|. This 
distortion, as well as the small orthorhombic deviation 
from the cubic symmetry, were not even observed in Refs. 
[ToL ITU Since these distortions are small, several theoret- 
ical papers chose to ignore them and assume cubic sym- 
metry. The corresponding cubic model has been taken 
as the "minimal" model needed to explain the physics 
in these materials. In the present Letter we show that 
although this model is of great theoretical interest, it is 
insufficient to explain the experiments. 

For the cubic titanates, there is one d electron in the 
t>2g degenerate manifold, which contains the wavefunc- 
tions \X) = d yz , \Y) = d xz , and \Z) = d xy . The large 
degeneracy of the resulting ground states, which involve 
both the spin and the orbital degrees of freedom [3L lil Il2^| . 
may then yield rich phase diagrams, with exotic types of 
order, involving a strong interplay between the spin and 



orbital sectors (e.g. 13. [ToL TTTh . justifying the broad the- 
oretical interest in this cubic limit. As we show, the cor- 
responding superexchange Hamiltonian (hereafter called 
the cubic Kugel-Khomskii (KK) model|T3) contains sev- 
eral interesting hidden symmetries. In addition, our anal- 
ysis shows that the KK Hamiltonian cannot yield some 
of the predictions which were claimed in the literature 
to follow from it. In particular, it has been suggested 
[l3| that the KK Hamiltonian gives rise to an ordered 
isotropic spin phase at non-zero temperatures, and that 
an energy gap in the spin excitations can be caused by 
spin-orbit interactions . We use the symmetries of the 
KK Hamiltonian to show that both of these predictions 
cannot hold. The observed long range order and finite 
gap [lCj must therefore be based on more complicated 
Hamiltonians, which go beyond the scope of this paper. 

As can be seen from Fig. ^ cubic symmetry implies 
no hopping (via oxygen p-states) among orbitals of type 
I a) along the a-axis. Neglecting the direct Ti-Ti hop- 
ping, KK called this axis the "inactive axis" for a or- 
bitals. This statement forms the basis for the remarkable 
symmetry properties of the KK Hamiltonian reported in 
this Letter. Apart from constant terms, the perturbative 
expansion of the Hubbard Hamiltonian with hopping of 
this type, to order e, yields the cubic KK Hamiltonian, 
H = TL X + H v + Ti z , where 
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and (ij) € a denotes a nn bond along the a-axis. Here, 
creates an electron at site i in a (3 orbital with spin 



-i,p,(T 



a, and one assumes that there is exactly one electron on 



each site, i.e. J2 a '• 



1, with n ia = J2 



„t 



For some purposes, it is convenient to separate the spin 
and orbital degrees of freedom^]- Defining the spin of 



2 




FIG. 1: A schematic view of the \Z) = d xy orbitals and 
the (indirect) hopping parameter t via intermediate oxygen 
p-orbitals. Positive (negative) regions of wavefunctions are 
represented by dark (light) lobes. Symmetry forbids indirect 
hopping along the z-axis for an electron in the Z-orbital. 

an electron at site i as Si, one has 

W a = (e/2) £ :1 • IS. • S, :• ./''. , (2) 

(ij)ea 

Jij = n ix n jx + n iy n jy + a\bibtaj + //</,<•/ J'/, . (3) 

Here, a\ and b\ create spin-less electrons in orbitals \X) 
and \Y), respectively, and rii X — a\a,i etc. 

Both Eqs. and iJSJ imply that whenever an a- 
orbital electron is destroyed, an a-orbital electron is cre- 
ated on either the same or another site. Therefore, the 
total number of electrons in each orbital is a good quan- 
tum number: any eigenfunction can be labeled by the 
total number of electrons in each orbital (i.e. Nx,Ny, 
and iV?)|l5|. Furthermore, when an a-orbital electron 
is destroyed, it is replaced by another a-orbital electron 
in the same plane perpendicular to the inactive (a) axis. 
Thus, for the n'th plane perpendicular to the a-axis, the 
total number N na of electrons in orbital |a) is conserved, 
i. e. it is a good quantum number. For example, for a 
cube of eight sites, the numbers Nix and N 2 x, which re- 
spectively are the numbers of X-orbital electrons in each 
of the two planes perpendicular to the x-axis, are con- 
served, and similarly for y and z. Thus, the states of 
the cube can be labeled by the six quantum numbers, 
{N 1x ,N 2 x,N 1 y 7 N 2 y 7 N 1 z,N 2Z ). 

Remarkably, there are more conserved quantities, as- 
sociated with electron spins. Defining the spin of an elec- 
tron in orbital a at site i as S ia = J2 pn c l,a, P i^]pv c i,a^/ 2 ^ 



where er represents the vector of Pauli matrices, and the 
total spin of all such electrons - located in an arbitrarily 
chosen plane #n perpendicular to the inactive a axis - 
as iS raQ = J2i£n &ia, we next perform a uniform, but ar- 
bitrary, rotation of all the S !Ct 's with i£n: we introduce 
an arbitrary 2x2 unitary matrix U^ n \ and write 

4, a ,„ = £^<<*,n' ien. (4) 
n 

Electrons in other orbitals or in other lattice planes are 
not affected by this transformation. Substitution of this 
transformation into Eq. shows that it leaves H in- 
variant. As a consequence of this symmetry, if one as- 
sumes long-range spin order, the spins associated with a- 
orbitals within any given plane can be rotated at zero cost 
in energy, thereby destroying the supposed long-range or- 
der. It also follows that S na commutes with TL, and thus 
both |iS rlQ | 2 and [5 nQ ,] 2 are good quantum numbers for 
each value of n or a. These symmetries can also be ob- 
tained from the original Hubbard model, provided one 
neglects Coulomb exchange interactions. 

This situation allows a rigorous proof of the nonexis- 
tence of long-range spin order at any nonzero tempera- 
ture for 7i of Eq. . Following the procedure of M ermin 
and Wagner (MW) . we choose 

C = 5+(k)=£ e - k - R c ^ QT c R , Q , i , 

R 

A = 5-(k + K)=£ e i; ( k+K )-% Qa c R , Q , T , (5) 

R 

where K is the wavevector of the order we wish to dis- 
cuss. Here aK = (n, it, it) is most relevant. Assuming 
long range order of S^(K) = s a . z and a corresponding 
staggered field h (for spins in the a-orbital), we end up 
with the MW-like bound 

1 > 2fcTK z | 2 l£[/i SQ , 2 + J Q (k)]~\ (6) 

k 

where J Q (k) oc Ylgd a {^ ~~ e~ jk '' 5 ) is proportional to 
the k-dependent parts in the Fourier transform of the 
non-zero nn spin-spin interaction (with nn vector 5) in 
the Hamiltonian. Since spins in orbital a couple only 
within planes perpendicular to the a-axis, it follows that 
J*(k) cx 2X^11 - cos(fc /3 a)] « J2p^ a a2k = a 2 k 2 ± ? , 
with no dispersion in the a- direction. For systems in 
d < 3 dimensions, the sum in Eq. © diverges as h — > 0, 
implying that s QjZ must go to zero. The conclusion is that 
the KK model is at its lower critical dimension d< = 3 
and does not support long-range spin order at T > 0. As 
we show elsewhere, a similar proof can be formulated for 
the original Hubbard model p7j. 

The same conclusion also follows from a renormaliza- 
tion group analysis of the model at finite T[jj}. Gen- 
eralizing to m orbitals and n-component spins, the spin 
free energy functional maps onto that of the "canonical" 
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FIG. 2: Spin and orbital configurations for a cube of eight 
sites. The thick lines indicate singlet spin states (dimers) and 
the X and Y indicate the orbital states of the electrons. Con- 
figurations in (a) and (b) are the dominant ones in the ground 
state wavefunction. The less dominant configuration (c) is ob- 
tained from (b) by allowing the interchange of two (X and Y) 
electrons along the z-axis, retaining their membership in the 
spin singlets (even though their positions have changed). 

nrn-component spin problem |18| . but with a (d — 1)- 
dimensional transverse gradient term, 

F = ^E( r + «!«)S a (q)-S a (-q) 

qa 

R ^ a a</3 ' 

where S Q (q) is the Fourier transform of S Q (R). Similar 
forms arise in connection with Lifshitz-like behavior |19j. 
This anisotropic gradient term shifts both the upper and 
the lower critical dimensions up by 1. For 3 < d < 5 
dimensions and n > 1 it yields decoupled n-component 
critical behavior. The free energy Q also reflects the 
symmetry with respect to independent rotations of the 
spin S Q associated with the single orbital a. 

These symmetries and conservation laws are very use- 
ful in the exact numerical diagonalization of finite Ti clus- 
ters, which indeed confirms their validity. We demon- 
strate this for a cube of 8 sites (Fig. |2J). Since the 
Hamiltonian commutes with the total spin S = £\ Si, 
the eigenstates can be identified by the quantum num- 
bers S [where S 2 = S(S + 1)] and S z . Since the energy 
does not depend on S z , it suffices to study the subspace of 
70 x 3 8 = 459270 states with S z — 0. A numerical analy- 
sis of the low-energy spectrum of this huge sparse matrix 
yielded 3 degenerate S = ground states, ^f y , and 
^> z , related by cyclic permutations. ^f z has N x = N y = 4. 
This information suffices to find the manifold containing 
each of these ground states. Since ^> z , for instance, is 
not degenerate within the manifold N x — N y — 4, it 
must also have the quantum numbers iVijc = N2X = 2, 
Niy — N2Y — 2, and N\z — N2Z — 0. A nonsymmetric 
choice like Nix — 3 and N 2 x = 1 would be degener- 
ate with Nix = 1, N2X = 3. The lack of degeneracy also 
implies that the total spin of the N na = 2 electrons in or- 
bital a in the n'th plane perpendicular to the a-axis must 
be S na = 0. Examples of such configurations, containing 
dimers of a-orbital electrons in the ev-planes, are shown 
in Fig. |21 The Hamiltonian allows an exchange of an X- 



electron with a Y-electron only along the z-axis (the only 
axis along which both types can hop) . Starting from Fig. 
Ufa), and performing all such possible exchanges, creates 
a manifold of 16 states (3 of which are shown in the fig- 
ure). Two other states with the same N na 's, but with 
the dimers along the z-axis, form another manifold, of 
higher energy. Indeed, a diagonalization of the resulting 
16 x 16 matrix reproduced the same ground state energy 
as found from the 459270 x 459270 matrix, demonstrating 
the power and the correctness of these symmetries. We 
are currently extending these numerical studies to even 
larger systems (such as N = 16 sites) to better under- 
stand the nature of the ground state in a real system. 

Since the KK Hamiltonian (|TJ forbids long range spin 
order at T > 0, the existence of such order (as in 
LaTiOs |10|) must result from some additional mecha- 
nism. Even for cubic symmetry, such mechanisms could 
include the small direct Ti-Ti hopping along the inac- 
tive axis, Coulomb exchange terms in the original Hub- 
bard model, or the spin-orbit interaction. In the real 
orthorhombic titanates one must also include JT distor- 
tions and oxygen octahedra rotations. A full discussion 
of all these effects lies beyond the present paper. Since 
the present paper concerns mainly symmetry arguments, 
we concentrate here on adding the spin-orbit interaction 
to the cubic KK Hamiltonian, where we can use such ar- 
guments to show the absence of a spin gap. Specifically, 
the spin-orbit Hamiltonian is 

iaf3aa' 

where L^g = (a 1 1/^1/3) is the orbital angular momentum 
matrix element. Since Tiso mixes orbitals, an cv-electron 
can hop via Tiso to orbital j3, then hop to a /3-orbital on 
a nn along the a-axis, and finally use TLso to return to 
orbital a. This generates an effective hopping between 
cv-orbitals along the "inactive" a-axis, invalidating the 
above arguments, shifting the lower critical dimension 
for total spin ordering back to d < = 2 and restoring long 
range spin order at d = 3. This mixing also eliminates 
the independent symmetries which we found for electrons 
within each orbital separately. However, as discussed be- 
low, there still remain some global symmetries for the 
total spin. Based on the signs of the leading couplings, 
we assume that the total spin orders antiferromagncti- 
cally, and proceed to show that the spin wave excitations 
in the ordered phase must be gapless. 

Again, an exact symmetry analysis clarifies the situa- 
tion. For electrons within each of the three degenerate 
t 2g orbitals discussed here, we introduce the following 
canonical transformation^ from spin to pseudo-spin: 

c* = VU 1 " 1 / (9) 

n 

where U a — <r a represents a different rotation for spins 
in different orbitals. As discussed in Ref. |6j, all the terms 
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in both the original 7i and in TLso now contain only com- 
bined operators of the form X^<4a ^,0, m> coeffi- 
cients which do not depend on the pseudo-spin indices /i 
(see Eq. (6) in Ref. |6J). These combined operators, and 
therefore also the full Hamiltonian, are rotationally in- 
variant in pseudo-spin space. Said differently, the Hamil- 
tonian is invariant with respect to a transformation on 
the original spins of the form c f wi = £„, V^JcJ^, 
with V' Q ) = ir Q U(T a , where U is an arbitrary unitary 
matrix. Thus the system possesses a continuous sym- 
metry, but it is not the usual symmetry with respect to 
rotation of the total spin. In the antiferromagnetically 
ordered phase, the spin staggered moment selects an ori- 
entation, and therefore the pseudo-spin will also exhibit 
broken symmetry. Rotation of the pseudo-spin will give 
rise to a manifold of zero energy states. This continuous 
symmetry guarantees that we have a (probably propagat- 
ing) zero-energy hydrodynamic mode : 20j. The rigorous 
conclusion is then that spin-orbit interaction permits the 
existence of long-range order at nonzero temperatures, 
but does not cause a gap in the elementary excitation 
spectrum, contrary to the assertion in Ref. Il4l Since our 
argument is based on symmetry considerations, it holds 



no matter what type of fluctuation is considered, and re- 
gardless of the orbital ordering (long ranged or liquid) . In 
analogy with results of Refs. Irll7l l2ll it is probable that 
when Coulomb exchange interactions and/or canting of 
the Ti-O-Ti bonds are included, spin-orbit interactions 
would lead to an energy gap in the excitation spectrum. 

In conclusion, we uncovered several novel symmetries 
of the KK Hamiltonian for cubic ti g systems. It is sur- 
prising that the KK Hamiltonian has been widely used 
in the study of interesting spin-orbital physics of transi- 
tion metal oxides for a long time but yet its remarkable 
symmetry properties were missed until now. Using these 
symmetries, we rigorously showed that the KK Hamilto- 
nian without spin-orbit interactions does not permit the 
development of long-range spin order in a three dimen- 
sional cubic lattice at nonzero temperature. Inclusion of 
spin-orbit interactions allows the formation of long-range 
spin order, but the excitation spectrum is gapless. 
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